Perturbative nonequilibrium dynamics of phase transitions in an expanding universe 
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A complete set of Feynman rules is derived, which permits a perturbative description of the nonequi- 
librium dynamics of a symmetry-breaking phase transition in \(f> 4 theory in an expanding universe. 
In contrast to a naive expansion in powers of the coupling constant, this approximation scheme 
provides for (a) a description of the nonequilibrium state in terms of its own finite- width quasipar- 
ticle excitations, thus correctly incorporating dissipative effects in low-order calculations, and (b) 
the emergence from a symmetric initial state of a final state exhibiting the properties of sponta- 
neous symmetry breaking, while maintaining the constraint (0) = 0. Earlier work on dissipative 
perturbation theory and spontaneous symmetry breaking in Minkowski spacetime is reviewed. The 
central problem addressed is the construction of a perturbative approximation scheme which treats 
the initial symmetric state in terms of the field <j), while the state that emerges at later times is 
treated in terms of a field f, linearly related to <f> 2 . The connection between early and late times 
involves an infinite sequence of composite propagators. Explicit one-loop calculations are given of 
the gap equations that determine quasiparticle masses and of the equation of motion for ((f) 2 (t)) and 
the renormalization of these equations is described. The perturbation series needed to describe the 
symmetric and broken-symmetry states are not equivalent, and this leads to ambiguities intrinsic to 
any perturbative approach. These ambiguities are discussed in detail and a systematic procedure 
for matching the two approximations is described. 
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I. INTRODUCTION 

It has long been apparent that the hot, dense mat- 
ter present in the early universe might undergo a variety 
of phase transitions [IJ. In particular, the "new infla- 
tion" scenario of Linde and Albrecht and Steinhardt 
j| (modifying an earlier proposal of Guth suggests 
that a symmetry-breaking phase transition at the GUT 
scale would lead to a period of quasi-exponential expan- 
sion, perhaps with highly desirable cosmological conse- 
quences (see e.g. ||). Variants of this proposal are still 
of great current interest. Amongst cosmologists, the con- 
ventional view is that the state of a quantum field during 
an inflationary era is adequately described by classical 
field theory (except for the purpose of estimating den- 
sity perturbations that are taken to originate from small 
quantum fluctuations) and, on the basis of classical cal- 
culations, that the new inflation scenario does not work 
in detail. However, a fully quantum-field-theoretic anal- 
ysis of the dynamics of phase transitions in, say, grand 
unified theories, has never been given. 

Apart from its intrinsic theoretical interest, a solution 
of this problem is of some importance. If there is a cos- 
mological era in which matter is adequately described 
by a spontaneously broken gauge theory, then the con- 
sequences of the associated phase transition need to be 
correctly understood, whether they lead to the desired 
inflationary picture or not. Moreover, the above conclu- 
sions drawn from classical calculations are not necessarily 
secure. One reason is that the effective potential which 
traditionally appears in the classical equations of motion 
is not, in principle, an appropriate dynamical tool (for the 



reasons discussed in amongst others) . Another is that 
the value of a classical field (say, the expectation value 
of a quantum field) does not necessarily provide an ad- 
equate characterisation of the actual quantum state. In 
an influential paper §, Guth and Pi argue, for the case 
of an unstable free scalar field theory, that at sufficiently 
late times the quantum probability distribution evolves in 
an approximately classical manner, because the minimal 
quantum uncertainty is neglible. However, as these au- 
thors recognise, this does not by any means imply that, 
say, the expectation value of the energy density (p{4>)) 
that appears in the field equations of semiclassical general 
relativity is well approximated by p( (</>)). Furthermore, 
the era of quasiclassical evolution arises from a growth 
in the amplitudes of unstable field modes, which may be 
considerable (see also ||) so that the effects of interac- 
tions may become significant even in a weakly coupled 
theory. 

The purpose of this paper is to derive a set of Feynman 
rules, by means of which the nonequilibrium evolution of 
a quantum field theory may be estimated perturbatively 
during the course of a symmetry-breaking phase transi- 
tion. One would like, for example, to be able to solve the 
semiclassical field equations 



(1.1) 



where T^ v is the stress tensor of an appropriate quan- 
tum field theory, and the question we address is how the 
expectation value can be estimated. Here, we study the 
example of a self-interacting scalar field, but the meth- 
ods we develop can be generalised to the case where 
this field belongs to the Higgs sector of a spontaneously 
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broken gauge theory. For simplicity, we restrict atten- 
tion to a conformally coupled field in a spatially flat 
Robertson- Walker universe with the line element ds 2 = 
a 2 (t) \dt 2 — da; 2 ], so that t is conformal time and the 
spatial coordinates x are comoving. The action for this 
theory can be expressed as 



S((f>) = / dPx 



(1.2) 



where m 2 (t) = a 2 (t)rn^ and too is the bare mass of the 
corresponding Minkowski-space theory. With TOq > 0, 
we would expect the zero-temperature Minkowski-space 
theory to exhibit spontaneous symmetry breaking, and 
we take the initial state (at a time that we shall call 
t = 0) to be a state of thermal equilibrium at a temper- 
ature l//?o which is high enough for the symmetry to be 
unbroken. The expectation value of a Heisenberg-picture 
operator 0(t) is then 



(0(t)) = Tr \e- 0oH{o) O{t)\ /Ti 



(1.3) 



where H(t) is the Hamiltonian, which depends explic- 
itly on time through the time-dependent mass m(t). 
This choice of an initial state is somewhat artificial, and 
is motivated mainly by the fact that it yields a well- 
defined problem. It has, in particular, the property that 
(4>(x, 0)) = 0. This feature will cause considerable diffi- 
culty, since it implies that (4>(x, t)} = at all later times 
also, so it is worth discussing at the outset. Of course, 
this initial state is well-defined, and should evolve in a 
well-defined manner, which it is of theoretical interest to 
investigate. One may wonder, however, whether the con- 
dition ((/)) = is too special to warrant the technical ef- 
fort needed to deal with it, and we want to argue that it is 
not. We are supposing that, like the gauge symmetries of 
more realistic models, the symmetry (f> <-* —<f> is an exact 
symmetry of nature. This means that any field configu- 
ration 4>(x) is physically indistinguishable from the con- 
figuration — <p{x). There is therefore no physical mech- 
anism that can produce different probabilities for these 
two configurations, which implies that (4>(x,t)) vanishes 
identically at every spacetime point (x,t). (This does 
not by any means imply that the state is spatially uni- 
form, since a quantity such as the energy density, which 
is not constrained by symmetry, may perfectly well have 
a non- uniform expectation value even when (<fi) = 0.) 
In the context of semiclassical gravity, and in the spirit 
of the Copenhagen interpretation of quantum mechan- 
ics, it is reasonable to suppose that the density matrix 
can depend only on quantities which couple to (and can 
therefore be "measured" by) the classical spacetime. It 
might therefore depend on, say, the energy density and 
pressure, which might be spatially non-uniform, but not, 
independently of these, on ^, which cannot be measured 
either by the spacetime or by any other physical probe. 
If 4> belongs to the Higgs sector of a gauge theory, then 



these qualitative arguments are superfluous, since a well- 
known theorem of Elitzur Q assures us that its expecta- 
tion value must vanish. 

For a variety of reasons, some of which we shall have 
cause to discuss in detail later, perturbation theory is a 
rather limited tool for treating this problem, but we know 
of no other approximation scheme which might be used to 
treat the single scalar field considered here, or the more 
realistic models of particle physics that one might wish 
to tackle. The functional Schrodinger picture approach 
developed in jHj] seems to be restricted in practice to 
Gaussian wavefunctionals (and, perhaps, to scalar fields 
in planar universes). Motivated by the belief that grow- 
ing unstable modes make perturbation theory completely 
unreliable, Boyanovsky, Holman and de Vega, with sev- 
eral other collaborators, have studied in considerable de- 
tail the case of 0(A r )-symmetric scalar field theory in the 
limit N — ► oo (see jllj and references therein). In this 
limit, the problem can be solved exactly (up to the nu- 
merical solution of a set of integro-differential equations) 
without the aid of perturbation theory. While significant 
insights can be gained in this way, the large-iV limit is 
a rather special one (in which, for example, dissipative 
effects are absent). Corrections of order 1/N and beyond 
seem to be intractable, so this approach does not appear 
to provided the basis for a non-perturbative treatment of 
any more realistic models. 

In constructing a perturbative means of tackling the 
problem, we shall need to draw on earlier work which 
investigated firstly the possibility of treating dissipation 
by describing a nonequilibrium state in terms of its own 
quasiparticle excitations and, secondly, how the phe- 
nomena normally associated with spontaneous symmetry 
breaking can be recovered when (cf>) — 0, by dealing in- 
stead with ((f) 2 ) which is not constrained to vanish. The 
results of this earlier work are summarised in section || 
below. We find, in particular, that two different pertur- 
bative expansions are needed to describe the symmet- 
rical state which exists at early times and the broken- 
symmetry state which exists at later times. Section [II 
describes a construction of the path integral which facil- 
itates this dual description and derives the complete set 
of Feynman rules. In section IV, we apply these rules 
to determine, at the lowest nontrivial order of our ap- 
proximation scheme, the gap equations for quasiparti- 
cle masses and the evolution equation for (<fi 2 ). Both 
this equation of motion and the differential equations 
obeyed by the late-time propagators require boundary 
conditions, which we obtain from the continuity of ap- 
propriate expectation values. Renormalization of the gap 
equation, the equation of motion and the boundary con- 
ditions which apply to them is considered in section 
Finally, the virtues and shortcomings of the approxima- 
tion scheme we propose are discussed in section [vj]. 
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II. SUMMARY OF PREVIOUS RESULTS 

A. Dissipative perturbation theory 

Our calculations of nonequilibrium expectation val- 
ues are based on the closed-time-path formalism [jl2| [Tg) . 
More specifically, we adopt the path-integral technique 
described by Semenoff and Weiss il7j| , in which Green's 
functions are obtained from the generating functional 



z (ja) = / [d(j> a ] exp 



iS(4> a ) + i I d 4 xj 



(2.1) 



Here, the single quantum field <p(x,t) is represented by 
three path integration variables cf> a (x,t) (a = 1, • • • , 3), 
which can be envisaged as inhabiting three segments of a 
contour in the complex time plane. The segment labelled 
by a = 1 is associated with time-ordered products. It 
runs from an initial time (which we shall call t — 0) to 
tf —ie, where tf is the largest time in which we are inter- 
ested and e is infinitesimal. The segment a — 2, associ- 
ated with anti-time-ordered products, returns from tf — ie 
to — 2ie. Finally, the segment a — 3, which provides 
a path-integral representation of the initial equilibrium 
density operator at t = 0, runs from — 2ie to — i(3, 
where p is the inverse of the initial temperature. The 
action S appearing in Z(j a ) is 



d 6 x 



1 dt£(<t>x)- f * dt£((t) 2 ) 
o Jo 

+ i I dr £_e(0 3 ) 
Jo 



(2.2) 



where £(4>) is the orig inal Lagrangian density (in our 
case, that given in ( |l.2] )), while £e is the Euclidean ver- 
sion associated with the density operator, namely 



£ E (<h) = \ [dAz? + \ (v^ 3 ) 2 



(2.3) 



The source term in (2.1) is 



d x j ■ (f> - 



d i x{ / dt \j x {x,t)4> x {x,t)+j 2 {x,t)(j) 2 {x,t)] 

drj 3 (x,T)(f> 3 {x,T) } 



(2.4) 

Of particular importance are the real-time 2-point func- 
tions (a, (3—1,2) given by 



G a p{x,x') 







dja(x) djp(x' 

Tr[pT(4>(x)4>(x>))\ 
Ti[p$(x)$(x')] 



3a=0 

Tt[ P 4>(x')4>(x)] 
Tr[pf(4>(x)4>(x'))] 

(2.5) 



where T and T denote respectively time-ordered and 
anti-time-ordered products of the quantum field operator 
4>{x) and p is the initial density operator. Other expecta- 
tion values can, of course, be obtained from appropriate 
derivatives of Z(j a ). 

To evaluate these expectation values perturbatively, 
one splits the action into an unperturbed part So(0 a ), 
which is quadratic in (f> a , and an interaction Si n t(4> a )- 



S(4> a ) = So(4>a) + Sint(4>a) ■ 



(2.6) 



After an integration by parts, Sq can be written in terms 
of a differential operator T> a b as 

So&a) = -\J d4x 4>a{x)V ab {x, dMx) . (2.7) 

The lowest-order approximations to the 2-point functions 
are then the propagators g a b(x,x'), which are solutions 
of 



V ab (x,d^)g bc (x,x') 



g a b{x,x')T> bc (x' ,- d M /) 
-iS ac S(x - x') , (2.8) 



subject to appropriate boundary conditions, which arise 
from continuity of the </> a and their time derivatives 
around the time path, including the periodicity condi- 
tion (f>3(x,f3) = (f>i(x,0). The perturbation series for 
an expectation value of some product of fields can be 
represented in terms of the usual Feynman diagrams, in 
which propagator lines representing g ab connect vertices 
arising from S^t. Since we consider only a spatially ho- 
mogeneous system, we shall normally deal with Fourier 
transformed Green's functions 

G ab (t,t';k)= [ d 3 (x-x')e- ik < x -^G ab (x,x') (2.9) 



and with propagators g a b(t, t'; k) defined in the same way. 

The standard choice for Sq is simply the quadratic part 
of S. With this choice (and supposing, temporarily, that 
m 2 (t) < 0, so that there is no symmetry breaking), one 
finds that the g ab (t,t';k) are composed of mode func- 
tions fk(t), which are essentially single-particle wave- 
functions, together with constants which can be inter- 
preted as the occupation numbers of the correspond- 
ing single-particle modes. Because this perturbation the- 
ory has as its lowest-order approximation the theory of a 
gas of free particles, which do not scatter, the occupation 
numbers remain fixed at their initial values^]. In the prob- 
lem at hand, we deal with an interacting system, driven 



Tn an expanding universe, even a free quantum field theory 
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away from equilibrium by a time-dependent Hamiltonian 
and, under these circumstances, one would expect the 
state of the system to evolve in response to its changing 
environment. In standard perturbation theory, this state 
is represented at lowest order by the fixed occupation 
numbers. One therefore suspects that low-order calcu- 
lations of time-dependent expectation values (which are 
all that one can realistically hope to obtain) will be ade- 
quate only over a period of time which is short compared 
with some characteristic relaxation time. The relaxation 
effects which would cause occupation numbers (or some 
appropriate generalization of these) to evolve with time 
in the expected manner arise from the absorptive parts of 
higher-order loop diagrams, so the perturbation expan- 
sion is improved if these contributions can be partially 
resummed so as to appear in the propagators g a b(t, t';k). 

As explained in detail in fl20|-|22]| , this resummation can 
be achieved by a somewhat more sophisticated choice of 
the real-time part of Sq , namely 



V, 



So (01, 02) = So 

' ' dt 



(2), 



h,<h) 
d 3 k 
(2^ 



<f> a (k, t)M a p(k, t, d t )<f>p(-k, t) , 

(2.10) 



where is the quadratic part of S, <j) a (k, t) is the spa- 
tial Fourier transform of 4> a {x) and M. a p is a differen- 
tial operator chosen in the following way. The propaga- 
tors g a p(t,t';k) now obey the spatial Fourier transform 
of (2.8) with a differential operator T> a p that contains 
M. a p, and we would like them to mimic the full 2-point 
functions as nearly as possible. From the hermiticity of 
<j) and of the density operator p, it is straightforward to 
show that all the real-time 2-point functions can be writ- 
ten in terms of a single complex function H(x,x') as 

G a p(x, x') = Hp(x, x')9{t - if) + H a (x', x)9(t' - t) , 

(2.11) 

where H 1 (x,x') = H(x,x') and H 2 (x,x') = H*(x,x'). 
Clearly, we want g a p to have the sa me f orm, and the 
most general choice of T> a p for which ( |2.8| ) admits such 
a solution is 



reputedly exhibits a phenomenon described as particle cre- 
ation (see e.g. |ilj,|is|). This arises if, for example, one adopts 
a family, parametrized by a time t, of mode expansions of 
the field, with mode functions f^{t) which behave approxi- 
mately as exp[— iu>k(t — t)] when t is near i. Each choice of a 
set of modes in general requires a different set of occupation 
numbers nj^' to describe the same physical state, and these 
may correspond roughly to the numbers of particles detected 
by a comoving observer at time i. However, for a given choice 
of i, njp is fixed, and does not evolve with time t. 



aP 



3?H 
-"fkd t 



■0k- 



ia k 
+ ia k 



- lOL k 

ioik 



(2.12) 



where a k (t), k (t) and 7fc(i) are undetermined, real func- 
tions. These functions can be determined in a self- 
consistent manner through an appropriate renormaliza- 
tion prescription. Thus, with the choice ( |2.10|) for So, 
the interaction S- m t contains a counterterm —5 / <pAi(j> 
which we require to cancel some part of the higher-order 
contributions to G a p, thereby optimising g a p as an ap- 
proximation to G a p. In this work we will, in particular, 
use a renormalization prescription such that 0k(t) has 
the form fik(t) = fc 2 + M 2 (t), corresponding simply to a 
mass renormalization, though more general prescriptions 
are possible in principle. It is perhaps worth emphasis- 
ing that the complete action S is, of course, independent 
of M a p- The introduction of this counterterm does not 
change the full theory, but serves to optimise our choice 
of a lowest-order approximation. The structure of M a p 
is analogous to that of the effective action obtained, for 
example, in p3| by integrating out extra environmen- 
tal degrees of freedom. Here, nothing is integrated out, 
but we might say that each field mode is treated self- 
consistently as interacting with an environment provided 
by the remaining mod es. 

With T> a p given by J2.12 ), the solution to (2.8) can be 
written as 

g a p(t, t'; k) = hp(t, t'- k)9(t - t') + h a (t', i; k)9(t' - t) , 

(2.13) 

where hi(t,tf;k) = h{t,t';k) and h 2 (t,t';k) = h*(t,t';k), 
with 

h(t,t';k) = lexp(-±f*dt"j k (t")) 

x l(N k (t') + 1) f k (t)r k (t') + (N* k (t>) - 1) f k (t)f k (t')} . 

(2.14) 

The mode function f k (t) is a complex solution of 

+ /*(*)= 0, (2.15) 

satisfying the Wronskian condition 

/*(*)/**(*) -/*(*)/**(«)=-*. (2-16) 
while N k (t) obeys 



d t + 2in k (t) 



n k (t) 



•7fc(*) 



[dt +lk(t)]N k (t) 

= 2ia k (t), (2.17) 



where Qfe(t) = 1/2/fe (£)/£(£). This function is also re- 
quired to satify the subsidiary condition 



N k (t)+N* k (t) 



•2*n fc (t) [N k (t) - NZ(t)] 
•7fe(*)[JV*(t)+JV*(*)]=0, (2.18) 
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in order that dt [hk(t,t') — hk{t', i)]\ t=t , — —i, which re- 
flects the canonic al com mutation relation. This condition 
is preserved by fl2.17D , so it can be regarded as a con- 
straint on the initital values of Nk and Nk- These initial 
values are determined by continuity conditions around 
the closed time path, as discussed in detail in |Q . Qual- 
itatively, we see that 7)fc (t) is a damping rate for unstable 
quasiparticle excitations and that Nk(t) has a loose in- 
terpretation in terms of time-dependent occupation num- 
bers for quasiparticle modes. Indeed, with approxima- 
tions appropriate to very slow evolution, equation (2.17) 
reduces to a Boltzmann equation |pl| in which jk (t) and 
OLkit) provide the standard scattering integral. In this 
work we consider only a real scalar field, but the dissipa- 
tive formalism outlined here can be extended to complex 
scalar fields |24| and (with somewhat greater difficulty) to 
spin-i fermions | p5[ . The dressed propa gators obtained 
by somewhat different methods in p6|]27| ] to describe dis- 
sipation in systems close to equilibrium (or, at least, to 
a steady state) seem to be a special case of those given 
here. 

For future use, we note that, associated with T> a $ is 
an operator 



d t 



lk{t) 



d a0 (t,k)=[ ut 



(2.19) 



which has the property 



g(t,t";k)d(t",k)g(t",t';k) 

( -ig(t,t';k) , t > t" > t' 

= 1 ig(t,t';k), t' > t" > t (2.20) 

I , otherwise . 



B. Spontaneously unbroken symmetry 

As indicated above, we wish to follow the progress 
of a phase transition starting from a high-temperature 
state in which the symmetry (f> <-> — (p is unbroken and 
so (4>(x)) — at each point in space. Since the dy- 
namics is governed by a Hamiltonian which respects this 
symmetry, it is inevitable that (4>{x)) = at all subse- 
quent times. We nevertheless expect to encounter states 
in which the phenomena conventionally associated with 
spontaneously broken symmetry (that is, with a non-zero 
value of (4>)) are realised, and will refer to such states as 
having "spontaneously unbroken symmetry" . Here, we 
review the means by which a vacuum state of this kind 
in Minkowski spacetime can be treated |28| . How such 
a state may be seen to emerge from a phase transition 
is the central problem to be addressed in this paper and 
will be discussed in detail later on. 

Heuristically, we can envisage two candidate vacua, 
say |+) in which {(f>) = +cr/vA and |— ) in which (<fi) — 
— cr/vA- The overlap (+|— ) of these states vanishes ex- 
ponentially as the volume of space becomes infinite, so 



) also vanishes if O is any local operator or the 
integral over all space of a local operator. Consequently, 
if we consider a superposition \a) = y/ct\+) + \/l — a\— ), 
for which (a\(j)\a) = (2a - l)a/VX, then (a\0\a) = 
a(+\0\+) + (1 — at)(— \0\ — ). If the symmetry is exact in 
Nature, then any experimental probe can couple only to 
a symmetrical operator, for which (+|0|+) = (— \0\— ). 
In that case, (a|C|a) is independent of a, and no ex- 
periment can determine the value of (<fi). In statistical 
mechanics, we need incoherent sums over states such as 
|+) and |— ), but the same principle applies. 

At the formal level, we need a means of calculating 
expectation values of symmetrical operators without as- 
suming a non-zero expectation value for </>. In particular, 
we need a lowest-order Hamiltonian whose eigenstates 
yield the Fock space built on \a) rather than on |+) or 
|— ). To this end, we take the state of spontaneously un- 
broken symmetry to be characterised by the expectation 
value of <fi 2 , which is not constrained by symmetry. In- 
stead of the conventional field variable ip(x), defined by 



(j)(x) = o/\f\ + i/j(x) , 
we deal with a field ((x) defined by 
(f> 2 (x) = U 2 + 2UC(x), 



(2.21) 
(2.22) 



where U 2 = (4> 2 (x)). Taking m 2 (t) = TOq in (1.2), we 
find 

C= l -(l + 2U-\y 1 d,td»Q-±U 2 e 

+U [ml ^ ( + l -5\Q) In (1 + 2U^C) , 

(2.23) 

where the last term comes from the functional Jacobian 
of the transformation and an irrelevant constant has been 
dropped. Since (£) = 0, the linear term must van- 
ish to leading order, so we identify U = y/6m%/\(l + 
0(A)). The range of ( in the path integral is, of course, 
— U/2---00. But, since U is of order A -1 / 2 , the lower 
limit can be extended to —00 at the expense of correc- 
tions of order e" 1 ^, which do not c ontribute to pertur- 
bation theory. On expanding (2.23) in powers of A, we 
find 



1 



1 



£=-d^(--(2m 2 )( 



(2.24) 



At leading order, we see that particles created by £ from 
the spontaneously unbroken vacuum have the same mass 
\/2mo as those created from the broken-symmetry vac- 
uum by ip, and it is not hard to convince oneself 
that these two types of particle are completely indistin- 
guishable. In this formulation, interactions arise from the 
kinetic term, and there are an infinite number of vertices 
(though only a finite number of these appear at a given 
order in A) . Renormalizability follows from the fact that 
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4> 2 {x) is a multiplicatively renormalizable composite op- 
erator in the standard formulation, and the 8 (0) in the 
Jacobian serves to cancel quartic divergences arising from 
the derivative interactions. 

It will be of some importance in what follows that the 
above manipulations do not, in themselves, assume that 
(0) = 0, but rather leave (0) undetermined. We may 
therefore consider the case that (0) = cr/VA and ( is 
related to tp by 

1 



C(z) = [<r 2 M - U 2 + 2 (a/Vx) rl>(z) + tp 2 (x) 



(2.25) 



In particular, we can express the connected two-point 
function (((x)((y)) c as 

(C(x)C(»)>c= ^<*(^(y)}c + ^j- 2 ^{x)^ 2 {y)) c 

+ 27^T ^ 2{X)HV))C + 7772 ^ 2 W^ 2 (y)) c • 



4t/ 2 



(2.26) 



In the Minkowski-space theory, one can easily calculate 
U in terms of a and verify this relation order by order 
in A. In the noncquilibrium state resulting from a phase 
transition that we plan to investigate, the field tp is not 
defined. Nevertheless, it will be possible to identify by 
eye term s corr esponding to the various correlators on the 
right of (HH). 

In what follows, we shall consider only the theory of 
one real scalar field. It is appropriate to mention, how- 
ever, that the analysis of spontaneously unbroken sym- 
metry outlined here can be extended to the Higgs sector 
of a spontaneously broken gauge theory [^9|. In accor- 
dance with Elitzur's theorem M, non-zero expectation 
values need be assigned only to operators which are in- 
variant under the gauge and global symmetries of the 
theory. The generation of fermion masses can also be 
accomplished in this way. Consider, for example, a pair 
of massless fermions ipL and ipR and a complex scalar 
with Lagrangian density 



- f (i>Li>R<f> + 4>R*pL<P*) + ' 



(2.27) 



With (0*0) = U 2 , we write <j>(x) = y/U 2 + 2U((x)e i6 W 
and construct the single massive spinor 

^(x) = e^W^i) + e %e{x)/2 tP R {x) (2.28) 

to obtain 



C = dpCWC + tp(i$ - fU)i> + (U 2 + 2UC)d^ 
+ ^^SJ + • • • . 



(2.29) 



Since only the derivatives of 9 appear, we need not assign 
8 an expectation value, and the fermion mass is fy/ (0*0). 



This indicates both that our formalism can in principle 
be extended to more realistic particle-physics models and 
also that our earlier remarks concerning the unobserv- 
ability of (0) are not invalidated by the linear Yukawa 
coupling of to fermions. 



III. DERIVATION OF THE FEYNMAN RULES 

The central purpose of this paper is to derive a set of 
Feynman rules which will permit us to follow the evolu- 
tion of the state of our system as the symmetry becomes 
spontaneously unbroken, in the sense described in the 
last section. As an intuitive guide to the considerations 
involved, we offer in Figures la and lb an artist's impres- 
sion of a time-dependent effective potential which might 
be thought to govern the evolution of 0. At times earlier 
than, say, t\, it has a single minimum at = while 
at times later than t\ it has two symmetrically placed 
minima at = ±0o(£). This is no more than an intu- 
itive guide, because we have given no precise definition 
of V e g ((/>). There are several conventional definitions of 
effective potentials, all of which refer to equilibrium sit- 
uations, and do not necessarily fit the dynamical role 
which is often forced upon them. 

Figures lc and Id show an artist's impression of the 
probability density V(<p,t) for the field at the spatial 
point x to have the value at time t. Because of spatial 
homogeneity, this probability is independent of x. Again, 
this gives only an intuitive expectation. Perturbation 
theory gives no ready access to V(4>,t) and its structure 
might well be more complicated than that sketched in 
the figure (and, of course, we are ignoring the problems 
of renormalization that would be encountered in trying 
to construct a well-defined V((f>,t)). Also, V(<f),t) does 
not by any means give a complete characterization of the 
state of the system. For that, one would need the vastly 
more complicated probability density for field configu- 
rations over all space. The intuition illustrated is that 
before some time to i the most likely value of cj) is zero, 
whereas after to, a broken-symmetry state emerges in 
which the most likely values are ±<r(t)/VA. 

These heuristic considerations do not provide a formal 
basis for the calculational scheme we wish to propose. 
The foregoing discussion is intended to motivate the as- 
sumption we do make, namely that there is some time 
to before which the most likely field values are near zero, 
and the appropriate field variable for perturbative calcu- 
lations is 0, while after to the most likely values are near 
±<r(t)/y/\ and the appropriate field variable is £. For 
brevity, we shall refer to these two varieties of perturba- 
tion theory as the "0-theory" and the "(-theory" . This 
situation presents two difficulties. One is that in the 0- 
theory the perturbative expansion of (0 2 ) has a leading 
term of order A , while in the (-theory the leading term 
is of order A -1 . Thus, although both perturbation se- 
ries are expansions in powers of A, they are not the same 



G 



expansion. We will propose a concrete means of dealing 
with this later on. The second diffic ulty is that we need 
to make the change of variable ( 2.22| ) only for times later 
than to, and thus to evaluate a path integral of the form 



t<t 



[d(]e 



(3.1) 



t>t 



This does not factorise into two independent integrals, 
on account of the boundary condition that 4> 2 (to) = 
U 2 (to) + 2U(to)((to)- Indeed, the nonequilibrium state of 
£ at time to is determined by the evolution of <fi between 
t = and t = to, and must be incorporated through a 
correct handling of the boundary conditions at to- The 
evaluation of this path integral is a hazardous under- 
taking, and the method we propose is not entirely rig- 
orous. We shall encounter various ill-defined quantities, 
and sensible prescriptions for dealing with these will be 
needed. We therefore begin by explaining our strategy 
in the context of a toy model, for which we obtain what 
are manifestly the right answers. 



A. A toy path integral 

The toy model in question is simply a free scalar field 
theory in Minkowski spacetime. The real-time part of 
the action (including a source J Q (t)) is 



dt 



.h{t)4>i{t) + J 2 {t)4> 2 {t) 



(3.2) 



For notational clarity, we do not indicate explicitly the 
spatial integral or the gradient and mass terms, which 
play no direct role in this part of our analysis. Wc 
will evaluate the generating functional given by standard 
methods as 



Z(J a ) 



JS(<P,J) 



= const x exp 



■± / dtdt'J a (t)g a0 (t,t')J (t') 
Jo 



(3.3) 

by integrating first over 4> a (t) for t > to and then inde- 
pendently over the remaining fields. We write the sources 
as J a {t) — j a {t)9(to — t) + l a (t)9(t — to) and in the first 
instance set j a (t) = 0. The propagator g a p(t,t') now 
obeys ( |2.8| ) in the form e ai d 2 g 10 = —i8 a p5(t — t'), where 
eix = —£22 = 1 and ei2 = £21 = 0, and we have again 
suppressed the spatial derivative and mass terms. We 
define 



L a (t) 



r t f 



f dt'g af 3(t,t')lf3(t') 

J t 



(3.4) 



which clearly satisfies e a pd 2 Lp(t) 
the change of integration variable 



-il a (t), and make 



<f> a (t)^<f> a (t) + iL a (t)0(t-ta), 



(3.5) 



After an integration by parts using the boundary condi- 
tion 4>i(tf) = 4>2{tf), which also entails h{tf) = liitf), 
the action (|3.2|) becomes 



S 



(i 



e a 04>a(t)<^i3(t)dt + ie al3 L a (t)dt4>/3(t) 



t=to 



, / l a (t)g a p(t,t')lf3(t')dtdt' 

'to 

i[Lt(t )-Ll(t Q )]S(0), 



(3.6) 



provided that we set 9(0) = i. Here, the fields for t > to 
are decoupled from the sources, so we can integrate them 
out. To do this, we note that exp[iS(4>i, ^2)] arises from 
the product of two time evolution operators: 



[#](■ • .) e ^(^« = Tr [(• • ■)U^ 1 (0,t f )U l (0,t f )} 

= r Tr[(---)U- 1 (0,to)U^ 1 {to,t f )U l (to,t f )U 1 {Q 1 to)}, 

(3-7) 



where Ui is the time evolution operator in the presence 
of the source Ji , and that the derivation of the path in- 
tegral implies the boundary conditions </>i(i/) = 4>2(tf) 
and 4>i(tf) = 4>2(tf)- In the absence of sources between 
to and tj, the product U 2 ~ 1 (to,tf)Ui(to,tf) is the iden- 
tity. Consequently, the path integration for to < t <tt 
yields a factor of 1 together with the boundary conditions 
0i (^0) — 4>2(to) and (to) = ^2(^0) on the remaining in- 
tegral. 

Before evaluating the remaining integral, we reinstate a 
non-zero source j(t) for times before to- In the present ex- 
ample, this could have been retained throughout, at the 
expense only of additional terms in our equations which 
were irrelevant until now. When dealing with the dynam- 
ics of symmetry breaking, however, we will have more 
cogent reasons for introducing j(t) only at this point. 
The quantity still to be integrated is now exjp[iS'(4>, J)], 
where 



S' = 



2 e ap4'a 

(*)^(*) + (ia(*)+M*))*a(t) 



dt. 
(3.8) 



The quantity L a (t) is a distribution concentrated at t 
to such that for any pair of test functions f a (t) 



to 



L a (t)f a (t) = ie a pL a (t)d t fp(t) 



t=t 



(3.9) 



Evaluating the integral in the standard way, we obtain 
const x exp[— with 
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J = 



J a (t)+ L a (t) g a a(t,t') 



jp(t') + L f3 (t') 



dtdt' . 



(3.10) 



The propagator g a p(t,f) is a solution of the same equa- 
tion as g a p{t,t'). If we assume that these are the same 
solutions (which is true if they and their first deriv ative s 
coincide when t = to or t' = to), then we can use ( 2.2C| ) 
in the form 

g ai {t,t")e 1 sd t "gsf}{t" ,t) = —ig a p(t,i/) 

x [0(t - t")8(t" - t') - 6{t' - t")8(t" - t)} (3.11) 

to find 

dtdt'L a {t)g aP {t,t')j (t') 



dt [ ° dt'l a (t)g aP (t,t')j p (t') (3.12) 

to J0 







and 

rto 



dtdt'L a (t)g a p(t,t')Lp(t') 



= -i[L 2 1 {t )-L 2 (t )]S(0). (3.13) 



Combining these results with those in (3.6), we recover 
the expected result ( |3.3| ). In particular, the terms pro- 
portional to 5(0) cancel. It is straightforward to verify 
that these m anip ulations also work if we replace the free 
field theory (3^2) with the dissipative approximate the- 
ory fl2.10D , provided that the operator t a pd t in (3.6) and 



( |3.9D is replaced by d a p as defined in (2.19) 



B. Path integral for a symmetry-breaking phase 
transition 

We are finally ready to undertake our central piece 
of analysis, which is to derive a perturbative means of 
calculating the generating functional Z(j a , l a ) of Green's 
functions which involve 4> a (x,t) (a — 1,2,3) for < t < 
to and ( a (x,t) (a — 1,2) for t < t < tf. The derivation 
is quite lengthy, and we set it out in several steps. 

Step 1: The generating functional. In the first in- 
stance, we set j a — and define 



Z(0, l a ) = / [dcf>] exp 



iS(4> a ) 



+ i dt j dxl a (x,t)( a (x,t) 

J t 



, (3.14) 



where 



£ a (x,t) = [217C*)]- 1 [<t>U*,t)-U 2 (t)] 



and U 2 (t) — (cjy^ix, t)). This is the exact expecta- 
tion value which must, of course, be determined self- 
consistently by demanding that (£ a (x,t)) = 0. The 
reason for setting the source j a for 4> a equal to zero is 
that this source breaks the exact symmetry on which our 
manipulations depend. We shall introduce j a at a later 
stage, and subsequently discuss the exact meaning of the 
generating functional obtained by this route. 

Step 2: A change of variable. We now change the 
integration variables from 4> a (x, t) to Ca( x ^ t) on the whole 
closed time path, so as to avoid difficulties over the time 
derivatives of fields at t = to. The transformation is 
formally legitimate even for t < to, but for these times 
we shall eventually have to undo the transformation in 
order to construct a sensible perturbation theory. The 
resulting Lagrangian is rather complicated, owing to the 
fact that U(t) is now time dependent, and we give only 
the real-time part (depending on <j>\ and fo) of which we 
shall make explicit use. For later convenience, we write 
it as the sum of several terms: 



C — Co + Mj 



(i) 



C 



(>i) 



c 



tdct 



Ad + A; 



Of these, the first is 



(3.15) 



(3.16) 



where V a p has the form shown in (2.12), but with 0k(t) = 
k 2 + M 2 (t). Here and below, we deal with the spatial 
Fourier transforms of fields, propagators, etc., but will 
generally not indicate explicitly their dependence on k or 
the associated momentum integrals. The total derivative 
term, in which e\ = —€2 = 1, combines with the d 2 terms 
in —7i( a 'D a pC t [3 to produce ^< 
The second contribution 



In 1 



-c] - -c 



-UV aP + ( U + M 2 U ) e a p ( (3.17) 



contains the principal interactions. There are derivative 
inte ractions of the kind already encountered in section 
[II B| , and these have been expressed in terms of the oper- 
ator T> a p, so that use can be made of (2.8) in computing 



Feynman diagrams. There are also non-derivative inter- 
actions arising from the time dependence of U (t) . The 
linear countcrtcrm 



4i ) = - U -m 2 {t)U +{\/Q)U i e a ( a 



(3.18) 



appears because we have taken U 2 (t) to be the exact ex- 
pectation value of cf) 2 , which means that (( a ) = 0. This 
requirement will be implemented self-consistently by re- 
quiring the counterterm to cancel all higher-order contri- 
butions to (Ca)- Further counterterms, denoted by 
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c 



(>1) 



1 

2U 
1 

+ 4 



U + M 2 U- -U 3 

o 



(3.19) 



arise from the fact that we used a renormalized mass 
M(t) and the dissipative coefficients afc(t) and 7fc(t) in 
£o- Again, these quantities are to be determined self- 
consistently by using the counterterms to cancel appro- 
priate parts of higher-order contributions to the 2-point 
functions. The dissipative counterterm 



7<9 t + ^7 + ia 
-jdt — 57 + ia —ia 



-ia 
l , 



(3.20) 



differs from the M a p of section II A only insofar as it 
excludes M{t), which is treated separately.The two addi- 
tional counterterms 

Adc = X a V a 0C/3 — (^aVapXp + — Y a d a f}(p , (3.21) 

where X a (t) and Y a (t, k) are as yet undetermined func- 
tions, are total time derivatives. They will be used to 
facilitate the handling of boundary terms arising from 
integrations by parts in the evaluation of certain Feyn- 
man diagrams. More terms of this kind might be needed 
for calculations at higher orders than we consider explic- 
itly in this paper. The sum of terms given so far differs 
from the original Lagrangian by a total time derivative, 
namely 



C 



d 



td 



dt [4 

+ t a U Cq — X a d a p(,fj — Y a d a /3(p 



(3.22) 



At this point, £td could be integrated round the whole 
time contour to yield zero, when use is made of the 
boundary conditions at t = and t = tf. We do not 
do this, because we want to insert an extra boundary at 
to and to indicate explicitly how we treat the boundary 
conditions there. Finally, the Jacobian of the transfor- 
mation is provided by 



£jac = 2 54 (°) 



lnU + 77-Cl 



In ( 1 + 77(2 



(3.23) 



Step 3: The path integral for to < t < tf . Our next 
task is to evaluate the £ integral for times after to. As 
usual, it is necessary to extract the interactions and coun- 
terterms as derivatives with respect to the source. Thus, 
we write 

Z(0, 1 a ) = exp \iV ( {-iS/Sl a )\ Zx(l a ), (3.24) 
where — (C — Co) dt and 



Zi(L) = / K]exp 



rtf 

1 l Cdt +i (Co + l a Q a ) dt 



to 



(3.25) 



The notation J to dt indicates that the integral is over 
both the real-time segments of the time path for < t < 
to and the imaginary time segment. We shall assume that 
the contribution to from £ t d can be neglected. The 
integral of this contribution over the closed time path 
does indeed vanish, provided that the boundary condi- 
tions to be explained in Step 4 below are valid. We shall 
later argue on heuristic grounds that the total derivatives 
in jCtdc (which would also integrate to zero) ought nev- 
ertheless to be retained. Now the manipulations leading 
to (3.6) can be repeated, with the result 



Z(0,l a ) = exp [iV c (-i5/Sl a )] 



Here, is 



exp (-Ijf) / [d()e 



(3.26) 



J C = l a (t)gi (t,t')l p (t')dtdt' 

Jt 

+ie a pL a (to)Lp(t )6(0) . 



(3.27) 



where is the propagator for the field £, which we shall 
need to distinguish from that for 0, and 



M*) = / gi P (t,t')W)dt' . 

Jto 

The action in the remaining path integral is 

/*° /"*/ «-> 

Cdt + J C a dt + i L a (t)d a p(t)Cp(t) 



(3.28) 



t=t 

(3.29) 



provided, as w e assume, that any effect of the shift anal- 
ogous to (3.5) on the interaction terms in C can be ne- 
glected. There is, indeed, no effect if these interaction 
terms are taken to exist at all times less than, but not 
equal to to, but we are unable to prove that this is really 
legitimate. As in the toy calculation, the path integral 
for to < t < tf can now be performed, yielding a factor 
ofl. 

Step 4-' The path integral for t < to- We now trans- 
form the path integration variables from £ back to 0. In 
making this transformation, we assume that the bound- 
ary conditions Ci(*o) = £2(^0) an d Ci(*o) = Czfo) trans- 
late into 0i (to) = ^2 (to) and 0i(to) = 02 (to)- Indeed, we 
have already implicitly assumed that the same is true at 
tf. At the heuristic level, this seem justified if we regard 
the path integral as a sum over sufficiently smooth func- 
tions Q (t) or ( a (t) for which d4>\{t)/dt = 20 ct (t)0 ct (t). 
The action in the remaining path integral is just the stan- 
dard action for 0, except for the boundary term in ( |3.29 ). 







At this point, we add sources j a (t) for the fields 4> a (t), 
thereby defining 



[d<fi] e 



t<t 



■■ exp [iV C (-iS/Sl a )] 
x ex P (-i l 7 f ) 

dt^C +j a (t)<t> a (t) 

+iL a (t) [4> 2 a {t)~U 2 {t)}) 
where now L a (t) is a distribution such that 



Z (jai la) 

with 

S 2 = 



(3.30) 



Co 



L a (t)f a (t)dt=L fJ (t)d 0a (t) 



(3.31) 



(3.32) 



It should be apparent that Z(j a ,l a ) has the following 
significance: (i) Z(0,l a ) correctly generates the expecta- 
tion values we seek for t < t < tf, (ii) Z(j a , 0) correctly 
generates expectation values for < t < to, which we 
also need; (iii) however, because of the way in which j a 
has been introduced, Z(J a ,l a ) does not yield correctly 
the expectation values that would exist after to if a real, 
physical source jit) had been present at earlier times. 
This does not concern us, since for us j a (t) is merely a 
technical device for generating the expectation values of 
the source-free theory. 

Formally, the remaining path integral can be evaluated 
by standard methods, with the result 

Z(j a ,l a ) = exp \iV c (-iS/5l a ) + iV+(-i5/6j a j\ 

xZ (j a ,l a ), (3.33) 

where represents the interactions of cj) for t < to. 
Here, dissipative effects are taken into account by a coun- 



terterm 



2 Ya 



M 



in which M. a p has the same struc- 



ture as M. a a, but with different coefficients and 
7fc(i) associated with the behaviour of the <j) correla- 
tor. The functional Zo(j a ,la) is conveniently written as 
exp (— |/C) , with 



JC = JCq + K^c, + £tr + K-6 + ICo ■ 



(3.34) 



The various terms arise from different stages of the fore- 
going calculation as follows. The first, 



JCc = 



l a (t)g C al3 (t,t')l p (t')dtdt' . 



to 



(3.35) 



provides just the £ propagator as in (3.27). The second, 

ft 



K 



<K 



Ja(t)g* b (t,t';l) Jb (t')dtdt' 



(3.36) 



arises in the same way from the path integral over <p. 
However, the boundary term in (3.31), being quadratic 



in cf), leads to a modified differential operator T^tu(t;l), 
whose real-time components 



TrJt:l)=V 9 M)-i[ 1W " , 



Li(t) 



(3.37) 



depend on the source l a (t) for £(t). The propagator 
g^ b (t,t';l) is a solution of 



Vi c (t;l)gi b (t,t';l) = -iS ab S(t-t'), 



(3.38) 



and is therefore also a functional of l a (t) . For this reason 
too, the final Gaussian path integral that remains after 
the extraction of JCq and JC^ depends on l a (t) and on 
evaluating it we find 



Finally, 



Kt, = Tr In 



ICs = itapL a (to)L p(to)5(0) 



(3.39) 
(3.40) 



is the last term of ( 3.27 ) and 
/C = -]T I* dtL a (t)U 2 (t) 

= - [Lx(t) - L 2 (t)] (d t +^(t)) U(t) 



t=t 



(3.41) 



is the last term of (3.31). The superscript on ^(t) indi- 
cates that this is the damping rate for £ rather than for 



C. The Feynman rules 

In the usual way, perturbation theory n ow co nsists in 
expanding the interactions and in ( 3.33 ) in pow- 
ers of A, and we shall shortly describe the diagrammatic 
rules through which the terms of the perturbation series 
can be represented. First, however, we must confront a 
difficulty which was postponed earlier. Namely, the per- 
turbation theory for t < to , based on the field <f> yields an 
expectation value for (j> 2 which is of order A , whereas the 
scheme based on £ for t > to yields a leading term of order 
A -1 . Indeed, a perturbative treatment of the ^-theory is 
possible only if we take U(t) = yWW) = O (A" 1 / 2 ). 

We propose to resolve this conflict in the following way. 
For t > to, we define v(t) = A 1 / 2 ?7(t). Then the expec- 
tation value of a quantity A(£) can be estimated in the 
form 



(A(Q) = A™ [a {v) + Xax(v) + X 2 a 2 (v) + • ■ ■] , (3.42) 

which is a power series in A, provided that v is formally 
regarded as being of order A . For v(t), we will have an 
equation of motion roughly of the form 



v = fo(v)+Xh(v) + X 2 f 2 (v) + 



(3.43) 
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whose right hand side is again a formal power series in A. 
The desired expe ctatio n value will the n be obtained by 
substituting into (3.42) the solution of (3.43) which satis- 
fies the appropriate initial conditions. The initial values 
v(to) and v(to) are calculable in the ^-theory as power 
series in A, but are of order A 1 / 2 . In principle, each of 
these power series can be truncated at any desired order 
(or some partial resummation technique might perhaps 
be devised). However, the solution of (3.43) with the 
appropriate initial conditions cannot be expressed as a 
power series in A and nor, therefore, can the final result 
for (A(()). Nevertheless, this scheme leads to a sequence 
of approximations which in principle can be systemati- 
cally pursued to arbitrary orders. In view of the evident 
complexity of the noncquilibrium theory, though, a de- 
termination of the convergence or summability properties 
of this sequence is well beyond the analytical powers of 
the present author. 

With this approximation scheme in mind, we now de- 
scribe the Feynman rules for constructing the required 
power series. For t < t , the rules are just the stan- 
dard ones for X(f> 4 theory, with propagators constr ucted 



9abj(t,t' ,t") 



according to the prescription described in section [I A 



For t > to, things are more complicated: we discuss first 
the propagators and then the vertices. 



1. The propagators 

In standard perturbation theory, one has a lowest- 
order generating functional Zo(j a ,l a ) = exp[— i/C] in 
which K, is a quadratic functional of the sources, giving 
rise to propagator lines which connect vertices, and there 
is one such propagator for each particle spec ie s. He re, /C 
is the nonlinear functional defined by (3.34)-(3.41) and 
this gives rise to an additional infinite set of propagators 
which connect times before and after to. In effect, these 
propagators represent the nonequilibrium density matrix 
p(to) which, in the Schrodinger picture, would result from 
evolving the initial state from t = to t = to- 



Of the contributions to K listed in (3.34), the first, 
/C^, yields just the propagator g^ a p(t,t'), which is indi- 
cated by the solid line of Figure 2(a). The second con- 
tribution, JC^q involves g^ b (t,t'; I), which is a solution of 
(3.38). This solution can be written as 



n=0 J 

xg^(ti,t 2 )L yS (t 2 ) ■ ■■gt b (t n , t') , 

(3.44) 

where L a p is the array with L\\ = Li, L22 = L2 and 
L12 = L21 = 0. The first term of this series {n = 0) 
is just the </> propagator g^ b (t,t'), depicted by the bro- 
ken line of Figure 2(b). The term n = 1 is equal to 
J t / gaby(t, t', £")Z 7 (£")dt". The new propagator 



E 

a/3 



Jt(t,h)gt b (h,t') 



da/j(tl)4 7 (*l>*") 



ti=to 

(3.45) 



is depicted in Figure 2(c), where the open circle denotes 

f/(ii) _1 d(ti)\t 1= t ■ The terms n = 2 and n = 3 are 
shown in Figures 2(d) and 2(e). The dotted lines repre- 
sent g^(ti,tj) in which, after performing the derivatives 

in d, ti and tj are set equal to to- Explicit expressions for 
these propagators are straightforwardly obtained, but are 
somewhat cumbersome and will not be r eproduced here. 

In the same way, /Ctr defined in ( 3.39 ) is given (apart 
from an irrelevant constant) by 



E 



dt\ - ■ ■ dt n tr 



L{ti)g*(h,t2)L(t 2 ) 



x/(t 2! t 3 )i(*3)---5 (*n,*l)J , 

(3.46) 



The propagators corresponding to the first few terms 
of this series are depicted in Figures 2(g) - 2(i). For 
example, the first term is J t * g a (t)l a (t)dt, where 



g a {t) 



11(h) 



(X-7 C (ii)) 



gL(h,t) - g C 2a (h,t) 



(3.47) 



Here, use has been made of the fact that g^it^t) — 
h^(t,t ) is indepe nden t of a and /3, as is easily checked 
from ( |2.15 ) and ( ^.14 ). In this propagator, the two time 
arguments of g^ft^t') are set equal before acting with 

the derivatives in d. Later terms in the series for both 
K.^ and /Ctr involve derivatives of g^(t,t') whose time 
arguments are to be set equal after differentiation, and 
these equal-time limits are, unfortunately, not unambigu- 
ously defined. They can, of course, be made well-defined 
by specifying the order in which limits are to be taken, 
but we are not able to offer a well-motivated general pre- 
scription for how this should be done. We shall, however, 
describe examples of low-order calculations in which the 
appropriate prescription can be ascertained with reason- 
able plausibility. 



The contribution /Co to JC defined in (3.41) can be writ 
ten as f** g^(t)l a (t)dt, with 



9a(t) 



U(ti) (d tl -7%)) [gL(h,t) -g c 2a (tut) 



ti—to 

(3.48) 
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and is represented in Figure 2(f). It evidently has a struc- 
ture similar to that of ( |3.47 ), and the significance of this 



will become clear below. 

Compared with diagrammatic rules of the conventional 
kind, most of the propagators shown in Figure 2 look like 
vertices, so it is perhaps worth emphasising that they re- 
ally are propagators. They must be used to connect the 
vertices described below in order to form a valid diagram. 
For example, if the object in Figure 2(c) were a vertex, 
it could be combined with the propagator of Figure 2(b) 
to form a diagram with the topology of Figure 2(g), but 
this does not result in a valid expression. To make this 
apparent at the visual level, we envisage the propagators 
as "cables" terminated by "plugs" (indicated by the ar- 
rows in Figure 2) and will depict genuine vertices as pos- 
sessing "sockets" that accept such plugs. Then, a valid 
diagram can be wired up by plugging cables into sock- 
ets, but not by tying cables together. Formally, plugs 
corresponds to the sources j a or l a which accompany the 
propagators in Zo(j a , l a ), while sockets correspond to the 
derivatives in V^(—iS/Sj) and V*(— i6/Sl). Once a valid 
diagram is formed, account must also be taken of the 
spatial momenta which we have not indicated explicitly. 
These appear in the standard way, being conserved both 
at genuine vertices and at the circled vertices internal to 
the composite propagators. Each of these circled vertices 
contains the operator d(t; k) whose momentum k is that 
of the single ( propagator emerging from it. 



2. The vertices 



The interactions and in (3.33) are, of course, 
represented by vertices. Interactions of (f> exist only be- 
fore to while those of ( exist only after to and there are 
no vertices involving both fields. For <j>, we introduce 
a renormalized mass /j,(t) and a dissipative counterterm 



M-^git) as described in section II A. The interaction part 



of the Lagrangian density is then 



A_ 
4! 



(3.49) 



The vertices corresponding to the quartic interaction, the 
mass counterterm and the dissipative counterterm are 
shown in Figures 3(a), 3(f) and 3(g) r espec tiv ely. 

Interactions of £ are contained in (3.17) - (3.21) and 



( 3.23 ). In order to do perturbation theory, we set U(t) = 
\ ~ 1 /' 2 v(t) and expand in powers of A. For £i n t, we obtain 



A 



+ |A,- 3 (v + M\) ( Cl 4 -C 2 4 ) 
-fA^C 3 )^^--- • 



Cl) 



(3.50) 



The corresponding vertices are those shown in Figures 
3(b) - 3(e), where a stroke on one of the legs indicates 
the action of T>i a on the propagator attached to that leg. 



The linear counterterm (]3 .18 ) becomes 



= -A- 1 / 2 



m 2 v + ±v 3 ) 



^aCa 



(3.51) 



and is repres ented by Figure 3(h). The counterterms 
contained in (3.19) include a mass renormalization 



(i) + M 2 v ~ ±v 3 ) (Cf - (I) , (3.52) 



depicted in Figure 3(k), of which we shall make explicit 
use, together with counterterms associated with M^a- 
The first two of these are shown in Figures 3(1) and 3(p), 
but we shall not make e xplic it use of them. Finally, the 
two total derivatives in (3.21) are represented by Figures 
3(i) and 3(m), while the first two of th e inf inite sequence 
of terms representing the Jacobian in ( [3.23] ) are shown in 
Figures 3(j) and 3(n). 



IV. EQUATIONS OF MOTION AND 
CONTINUITY CONDITIONS AT to 

We now discuss several calculations which serve the 
dual purpose of completing the specification of the Feyn- 
man rules derived in the last section and of illustrating 
their use. The Feynman rules are so far incompletely 
specified for two reasons. One is that that they involve 
the expectation value v 2 (t) — A((j> 2 (t)}, the renormal- 
ized masses n(t) and M(t) and the dissipative coefficients 

a t (*)' a fc(^)> 7fc (*) an d 7fcM f° r which we have no con- 
crete expressions in hand. The second is that the prop- 
agators gf b (t, t'; k) and <?£g(t, f; k) are solutions of (2. J) 



(with the appropriate operator T> a b in each case), but the 
appropriate solutions must be identified through suitable 
boundary conditions. W hen t hese p ropag ators are rep- 
resented in the form of (2.13) and (2.14), the remain- 
ing ambiguity resides in the fu nctio ns Nf(t) and N^(t). 
These in turn are solutions of (|2.17) and it is the initial 



conditions for these functions that are needed. 

For the unbroken-symmetry state prior to to, the renor- 
malization and boundary conditions which determine 
fi(t), at(t), 7j£(t) and N%(t) are described in |^@] and 
we shall not repeat the discussion here. The determina- 
tion of aj,(t) and 7^(t) presents no new difficulty beyond 
that of computing the required integrals and will also 
not be discussed. It is therefore the determination of 
v(t), M(t) and the initial conditions on iV£(t) which are 
principally of interest. Also of concern is the fact that, in 
evaluating the path integral for Z(j a , l a ), we were unable 
to treat the boundary conditions at to in a fully rigorous 
manner, with the result that certain ambiguities remain. 
We will show how these ambiguities can be resolved at 
the lowest non-trivial order of our approximation scheme. 
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A. The condition (£(t)) = 



where 



We start with the fact that U 2 (t) = A~ 1 w 2 (t) was de- 
fined to be the exact expectation value of {4> 2 {t)) and is 
therefore to be determined self-consistently from the re- 
quirement that (C(t)) — 0. It is sufficient to ensure that 
the one-particle-irreducible contribution (£(£)) ipi van- 
ishes. As illustrated in Figure 4, this quantity can be 
separated into two parts, 



Cm - 



«;(*)> rpi = <c(*)>i a pT hored) 



<cw> (fr 



IPI 



(4.1) 



Diagrams contributing to the first part contain at least 
one of the circled vertices internal to one of the composite 
propagators, which we will describe as "anchoring" the 
corresponding time argument at t = to, whereas all the 
vertices in the second part are free to range between to 
and tf. In Figure 4(a), the diagram labelled (ii) is just 
the propagator of Figure 2(g), while that labelled (iii) is 
constructed from the propagator of Figure 2(c) together 
with the ordinary <f) vertex of Figure 3(a). These are the 
first two of a sequence of diagrams whose sum reproduces 
(</> 2 (io)) as calculated from the </>-theory. More precisely, 
on combining these with diagram (i) (which arises from 
Figure 2(f)), we obtain the expression 



(0 2 (*l)> 



(4.2) 



With the natural requirement that (4> 2 (t)) and its first 
derivative should be continuous at to, this vanishes iden- 
tically. The fact that contributions from several different 
terms of ( p. 34 ) conspire to give this satisfactory result is 
somewhat reassuring. We now turn to the second part of 
(CW)ipi illustrated in Figure 4(b). First, we dispose of 
the diagram labelled (v). In this diagram, acts on 
the internal propagator to produce S(0), and this contri- 
bution is precisely cancelled by the Jacobian counterterm 
of diagram (iii). In order to set (C(£))ipi equal to zero, 
we would like the remaining contributions to sum to an 
expression of the form J Zait^g^Jt' , t)dt' , so that z a (t) 
can be set to zero. The obstacle to this is diagram (vi), 
in which T>^ g acts on the external propagator to produce 
6(t — t'). Now, this problem could be solved through an 
integration by parts to make T) a o act on the bubble to its 
left, were it not for the boundary terms at t' — to, whose 
interpretation is a little unclear. It was to facilitate the 
handling of this integration by parts that we introduced 
the counterterm proportional to X a (t) in (3.21), which is 
a total derivative and contributes diagram (ii) of Figure 
4(b). In effect, we can now perform the integration by 
parts without incurring boundary terms by choosing 



Xl (t) = x 2 (t) = -^4(t), 



(4.3) 



d 3 k 



/i c (t,<;fc) 



(4.4) 



is (with the spatial momentum now made explicit) the 
bubble contained in diagrams (iv) and (vi). With this 
choice, the term Xa'D^g^p cancels diagram (iv) while 

—t^a'D'^gXp supplies the result of the integration by parts. 
The counterter m wa s, of course, subtracted off again in 
Ltd (equation ( 3.22| )), which we subsequently assumed 
could be neglected. This assumption is now seen to 
amount to a prescription for handling the present inte- 
gration by parts and other similar ones. In the present 
case, we take this prescription to be justified a posteriori 
by the fact that we can now set (((t)) = and hence ob- 
tain a sensible equation of motion for v(t). Including the 
remaining diagram (i) (which represents the counterterm 
Q3.18 )), the equation of motion at this order is 



v — m v H — v" 
6 




-± =0. (4.5) 



B. The 2-point function 

It is now necessary to address several issues con- 
cerning the connected 2-point function (i'))c- 
Clearly, if t < to and t' > to, this object is equal to 
(0a MC/3 (>')), while if t and t' are both greater 
than t , it is equal to 4U(t)U(t'){( a (t)(f)(t')). Continu- 
ity of this function at to will supply the initial conditions 
needed to specify (^(t, completely, and we also re- 
quire a suitable definition of the renormalized effective 
mass M(t). We shall deal with these issues at the lowest 
order of our approximation scheme, but first a technical 
question must be addressed. 



1. Anchored contributions to (C<*(£)C/3(0) 

The functional K, which gives rise to our v ariou s prop- 
agators contains a contribution ICs, given in ( 3.40| ) which 
is no t incl uded in Figure 2. In the toy calculation of sec- 
tion 



[II A, the analogous quantity was found to cancel 



exactly another singular contribution (equation (3.13)), 
provided that the propagators satisfied continuity con- 
ditions which were expected on other grounds. We now 
find a similar cancellation involving the composite prop- 
agator of Figure 2(h), which is the lowest-order anchored 
contribution to (( a (t)Cp(t')). The internal lines in this 
propagator correspond to the expression 
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d 3 k 



X%(t, t'; k) = J j^9%{t, t'; k')gl (t, t'; k' + k), R a (t; k) 



(4.6) 

and 2x^p{t,t';k) is precisely the lowest-order contribu- 
tion to (<^(i)0| (£')). Suppressing the momentum argu- 
ment, the expression represented by Figure 2(h) is 



gi-y {t,tx)d lS (*i)xt {h,t 2 )d e<7 (t 2 )gi f3 (t 2 , t') 



1— r 2— c 



As indicated above, this expression is ill-defined, because 



it involves the quantity c^c^Xafl^i, £2) 



How- 



■1— z 2— c O 



ever, we can use (3.11) (with e a pdt replaced by dap) to 
find that this contribution to /Ct r is exactly cancelled by 
ICs, provided that the boundary conditions on g^ a p{t,t') 
satisfy 



2X^(M') =0 AU{t)U{t')gi^t,t') 



(4.8) 



The notation =0 here indicates that the quantities on the 
left and right, together with their derivatives with respect 
to t and t' are equal at t — t' — to. If this prescription 
(which is consistent with the boundary conditions to be 
discussed below) is adopted, then we may delete the prop- 
agator of Figure 2(h). This matter will, however, need to 
be reconsidered when we come to discuss renormalization 
in section |v|. 



2. Continuity of the 2-point function 



At the lowest order of approximation, the 2-point func- 
tion (fiait)^^ 1 )) is represented by one of the diagrams 
shown in Figure 5, depending on whether its time argu- 
ments are greater or smaller than to- The anchored ver- 
tex in diagram 5(b) involves U(to) and U(to), which we 
expand in powers of A, retaining only the leading terms 



tf<°>(t)l =lf(t) 



d 3 k 



u = u(°Ht ); i>o= > (0) W 

dt 



(4.9) 



Because of the propagator structure exhibited in ( |2 . 1 3| ) 
this diagram can be expressed in the form 



gt%t';k) = -i 



d 3 k' 

(2tt) 3 



ht(h,t; k')ht(h,t; k' + k) 



(d tl - 7 c (*i)) [ft c (i',ii;fc)-ft c *(Oi;fc) 



*1— *0 

(4.10) 



In order to display the boundary conditions on the 2- 
point function succinctly, we define 



d 3 k' 



hi(t ,t;k')hi(t ,t;k' + k), 



(2tt) 3 
( d 3 k' 

S a (t;k)=d tl / ht{t u t-k')ht{t u t-k' + k) 



(2tt) 3 

P{t;k) = h c (t,t Q ;k) - h c * {t,t ;k) , 

Q(*;fc) = (fl ei -7^(*i)) 

x |V(£,ii;fc) - /i c *(t,ii;fc) 



tl=to 



(4.7) which satisfy, in particular, 



P(t ) = 0, Q(t )=i 
P(«o) = -i, Q(t )=0 

Ro(k) = fli(to;fc) = ik(*o;*) = R* {k) , 
S (k) = Si(t ;k) = S^(to;k). 

In terms of these, we have 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



U 



R a (t)Q(t') 



U 



+ \-^R a (t)-S a (t)\P(t') 



dt>(4>l{t)4>l(t>))\ t , =2S a (t). 



(4.15) 

(4.16) 
(4.17) 



In principle, the function {4> 2 a {t)4> 2 p{t')) should be per- 
fectly smooth as its arguments pass through to. We will 
actually attempt to impose the rather weaker conditions 
that the function itself and its first derivatives are con- 
tinuous at to. At lowest order, these conditions read 



&M) = ^rmt)^{to)), 



1 



2U 



gt$(to,t')=2U g^(t ,t') 



dtgi$(t,t>) 



t'=t 



t = t 



= 9v 



mt)m')) 



2l7(°)(t') 



2U^(t)gi p (t,t>) 



(4.18) 
(4.19) 



t'=t 
(4.20) 

(4.21) 



Of these, (|4.18| ) an d Q4.20D are satisfied i dentic ally, be- 
cause of ( [4.1 1|) and ( |4.12| ). The condition ( 4.19 ) is satis- 
fied provided that 



Im Sq 



-u 2 



(4.22) 



and 



NC (to) = 772 [^e So -{Z- ±7 C (M) Ro] , (4.23) 
^0 
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where 



(4.24) 



and /fc(t) is the mode function introduced in ( 2.15 ) that 
appears in g^. Explicit calculation shows that ( 4.22) ) 
is indeed true and (4.23) of course provides on e of the 
initial values that we seek. The last condition ( 4.2l|) is 
satisfied if 



^(*o) = 772 
^0 



Re S - 2ZReS 



Z 2 (t ) ) R 



with So = Si (to) = Sf(io), and if also 
ImS* = -f/ V(io). 



(4.25) 



(4.26) 



The second of our in itial values for the C prop agat or is, of 
course provided by ( 1.25| ), but the status of ( 4.26] ) is less 
clear. Evaluating this equation explicitly, we obtain a 
moderately plausible relation between the damping rates 
for (f> and £, namely 



d 3 k' 
(2nf 



7&(*o)**(to,*o; *:' + *) = U^i(to) . (4.27) 



However, this relation is not automatically satisfied. Ap- 
parently, it implies a constraint on the prescriptions used 
to define these damping rates, but we are not able to 
say whether a prescription of the kind described in pl| ] 
would naturally satisfy this constraint. In any case, when 
we come to consider renormalization, we shall find that 
the boundary conditions derived here require a significant 
modification. We note, though, that in this unrenormal- 
ized form, the initia l values N**(t Q ) and N**(to) do satisfy 
the condition (2.18). 



3. The 2-point function (( a (t)(p(t')) 

The last step required to specify the Feynman rules 
completely is to find a suitable definition of the renormal- 
ized C mass M(t). To do this, we calculate the 2-point 
function (Ca(t)Cp(t')) at first order in A. It is the sum of 
the eighteen diagrams shown in Figure 6 (recalling from 
our earlier discussion that the composite propagator of 
Figure 2(h) can be ignored). The computation of these 
diagrams is straightforward, except for Figure 6(s), in 
which T> a /3 acts twice on the same internal propagator. 
In order to evaluate this diagram we have (i) ignored 
some contributions from the dissipative coefficients a 1 * 
and 7^ which can formally be regarded as of higher order; 
(ii) adopted the following prescription for the equal-time 
limit of derivatives of the propagator: 



dtgL(t,t')\ 



ii m +ii m ) 9 *sL(M') 

Z I T t \,Z , 



\ [d t ht(t,t')+d t ,h C (t,t')]\ t=t ,; 



and (iii) identified the second counterterm in (3.21) as 



Y 1 (t)=Y 2 (t) = ~ 



Xlf(t) 



(4.28) 



where /£(t) was defined in (4.4), so as to facilitate an in- 
tegration by parts. The result we obtain is consistent, at 
thi s ord er of approximation, with the structure displayed 
in 

( g|| ). Specifically, we write it as 



1 



2 

A l/2 



1 



2 

A 1 



c 

v(t) al3 



(t,t 



4«(t) af} 



(t,f) 



v(t>) 

i 



v(t>) 



where the quantity 
a(t) - v(t) 



1 AJf(t) 
2 v 2 (t) 



0(A 2 



(4.29) 



(4.30) 



might loosely be interpreted as corresponding to the 
peaks of a probability density of the kind sketched in 
Figure 1(d). In the first term of (4.2E), with which we 
are principally concerned, we have 



T af3 



9U(t, t')+i df'gi^t, t")e lS A(t")g C 5(3 (t", f) 



to 

dt"dt'"gi y (t, t")e, 6 B Se (t", t"')e eX g{p(t"', t') , 

(4.31) 

with auxiliary functions A(t) and B a p(t,t') which will 
be specified shortly. It will be seen, for example, that 
the term (a/v) t g^ t p(t,t')(a/v)t' arises from the sum of 
diagrams (a), (g) and (h) of Figure 6, together with a 2- 
loop diagram which is easily seen to be present at order 
A 2 . As always in the closed-time-path formalis m, the 
signs e a j3 attached to the vertices and the structure (2.12) 
of the propagators ensure causality, in the sense that the 
integrands in (4.31) vanish whenever t" or t'" is greater 
than both t and t'. The auxiliary functions are given by 



A(t)- 



M 2 v 



| (^ ) ifw 



2v(t) 



m 2 ) [ h. 



9A 
2 



(4.32) 
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and 



B af3 (t,t')=-- 



9 A fv + M 2 v 



I2 (t, t') al 3 + ie a /35(t - t')h 



v + M 2 v 



(4.33) 



with 
rC/v +1. 



d 3 k' 



(4.34) 



The terms proportional to I2 cancel out in (4.31) and 
have been inserted for the purpose of renormalization, 
which is considered below. For completeness, we record 
that the remaining functions are 

G^f\t,t') = -3U 1 / 2 [* f dt"gUtX)e^ 



to 



G<*'*)((,t') = -3iA" 2 /' dt"4(t,t")a 



to 

v + M 2 v 

2 

V 



e -yS9s(}(t",t') 



G%* 2 \t,t')=2li{t,t'U. 

These are easily seen to c orres pond at lowest order to the 
correlators indicated in (2.26). 

Our renormalized mass M(t) will now be defined by 
requiring A(t) = 0. In the noncquilibrium theory, there 
is no clear analogue of the pole of a zero-temperature 
propagator in Minkowski spacetime, which unambigu- 
ously identifies the mass of a stable particle. Here, our 
rationale is, rather, to optimise g>£g as an approximation 
to the full 2-point function by using the mass countert- 
erm (3.52) to cancel as much as possible of the higher 



order corrections. Although the function A(t) was iden- 
tified by eye from our explicit O(X) result, we assume 
that this condition actually makes sense to all orders, so 
the definition of M(t) is an implicit one, which can be 
realised in practice only to a g iven order of approxima- 
tion. Setting the expression ( 4.32 ) to zero appears to 
give a second equation of motion for v(t) which issimi- 
lar, but not identical, to that already found in (|4.5| ) by 
requiring (((t)) = 0. In fact, these two equations to- 
gether yield both the equation of motion for v(t) and the 
relation between M(t) and the bare mass m(t). Consis- 
tently ignoring terms of order A 2 , we obtain the equation 
of motion 



(df + M 2 ) 



v-^l + 0(X 2 ) 
2 v 



y 



0(A 2 ) 



(4.35) 



and the gap equation 



m 2 (t) = -M 2 (t) + \v 2 (t) [l - AI 2 ] + 0(A 2 ) . (4.36) 

This gap equation expresses the bare mass m(t), which 
occurs only in the counterterm (3.51), as a function of 
M(t), v(t) and A. For a given m(t), these are two equa- 
tions to be solved simultaneously for v(t) and M(t). In 
the specific application to an expanding universe, we have 
m 2 {t) = a 2 (t)mp, and one must simultaneously solve Ein- 
stein's field equations for the scale factor a(t). 



V. RENORMALIZED EQUATIONS OF MOTION 
AND INITIAL CONDITIONS 

The equations of motion we have just derived in- 
volve integrals Ii(t) and /|(i, t'\ k) a p which are ultravi- 
olet divergent, and many more divergent integrals are of 
course to be expected at higher orders. Of necessity, the 
noncquilibrium theory is formulated in terms of quanti- 
ties which are defined implicitly as the solutions of dif- 
ferential equations or of self-consistency relations. This 
system of equations, while ultimately susceptible of nu- 
merical solution when truncated at some (no doubt low) 
order of perturbation theory, is far too complicated to 
admit of any exact analytical solution, and a proof that 
they can be renormalized so as to remove ultraviolet di- 
vergences at all orders is currently beyond the ingenuity 
of this author. A numerical investigation at the order 
of approximation we have considered explicitly ought to 
be feasible, provided that a suitable renormalization pre- 
scription can be given. Even at this level, we have in 
hand only formal expressions for the propagators g a b(t, t') 
in terms of mode functions fk(t) and generalized "oc- 
cupation numbers" represented by the functions (t) , 
whic h mu st be found from numerical solution of (2.15) 
and (2.17), so an unambiguous proof of renormalizabil- 
ity is rather difficult. Here, we suggest a renormalization 
prescription which might be applied in the context of 
a numerical computation involving a spatial momentum 
cutoff A, which reflects the renormalizations that are well 
known to work in the Minkowski-space theory. Since ul- 
traviolet divergences are a feature of the short-distance 
and short-time properties of the theory, one expects that 
they should be essentially independent of the nonequi- 
librium state, and that these renormalizations should be 
adequate. In practice, one hopes that numerical compu- 
tations would approach finite limits when A — > 00, and 
we shall offer circumstantial grounds for optimism that 
this should indeed be so. 
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A. Large-fc behaviour of propagators and 1-loop 
integrals 

We need to estimate the large-fc behaviour of our prop- 
agators, and do so using a method similar to that de- 
scribed, for example, in |30[ . The usual propagators of 
the equilibrium theory are obtained from our g a p(t,t'), 
by setting Nk(t) = 2nfc + l, where rik is the Bose-Einstein 
distribution, which falls off exponentially at large k. 
When time evolution is sufficiently slow, the dissipative 
coefficients ak(t) and Jk(t) are given by scattering inte- 
grals [^l] , which are convolutions of the n k , and also ex- 
ponentially small at large fc. For present purposes there- 
fore, we assume that it is sufficie nt to set Nk(t) — 1 and 
a k(t) = Jk(t) = 0. Then, from (2.14), we need only to 
estimate h(t,t';k) ~ hfk[t)fk{t')i with mode functions 
obeying 



M 2 (t)} f k (t) = 



(5.1) 



and the Wronskian condition ( 2.16) ) . The general solution 
may be written as 



/*(*) 



exp 



i I n k (t')dt' 
'o 



(5.2) 



where 9(k) is an arbitrary, but time-independent, func- 
tion of k which does not affect the propagator and 
will therefore be ignored. We will choose (as is always 
possible) a time-dependent frequency which behaves as 
= fc + O(fc ) for large k. In this case, we obtain the 
large-fc expansion 



2fc 



iw(t,ti)-^x(t,t>) + o(k-*) 



where 



W(t,t') 



M 2 {t")dt" 



(5.3) 



(5.4) 



and X(t,t') = M 2 {t) + M 2 (t') + \W 2 (t,t r ). This ex- 
pansion can be used to isolate the ultraviolet divergences 
of the two integrals l\(t) and I^t^t'^k). We impose a 
cutoff A on physical momenta fc p h = |fe p h| < A. Up to 
this point, our analysis has been entirely in terms of co- 
moving coordinates and momenta and, in principle, we 
should translate all of the above results into physical co- 
ordinates in order to impose this cutoff systematically. 
However, both if (t) and the divergent part of t'; k) 
are local in time and, for our present purposes, it is suffi- 
cient simply to convert the physical cutoff to a comoving 
cutoff fc < a(t)A. Specifically, we find 



i 

8^ 



(a(t)A) 2 - M 2 {t) In (a(t)A) 



(5.5) 



4 (t, t'; k) af} = -i5 aP e a S{t - t')I 2 (t) + 



— ln(a(t)A) 



(5.6) 
(5.7) 



where the ellipsis represents ultraviolet-finite contribu- 
tions . The quantity Ii(t') here is that appearing in ( 4.32| ) 
an d (|4.33|) and we see, in particular, that the net integral 
in (fh33|) is finite. 



B. Renormalized equation of motion and gap 
equation 

We hope, of course, that both the equation of motion 



(4.35) and the gap equation (4.36) can be expressed in 
a form which is free of ultraviolet divergences. To this 
end, we introduce the physical particle mass m, which 
locates the pole of the Minkowski-space propagator in 
the broken-symmetry vacuum. It is related to the bare 
mass mo by 



A 



16tt 2 



rA 2 



l 



rh 2 



1 



A 



16tt 2 



In 



A 



+ 0(A 2 ), (5. 



with c = 1 + In 2 — v3tt/2. (The 1-loop integrals used 
to obtain this relation are, of course, Minkowski-space 
integrals involving a physical 3-momentum fc p h, whose 
magnitude is cut off at the value A.) In the Minkowski- 
space theory, Green's functions involving the operator 
(ft 2 are rendered finite by a combination of additive and 
multiplicative renormalizations. Guided by the standard 
theory of these renormalizations, we surmise that a renor- 
malized expectation value v^(t) = (4 )2 )r can be defined, 
at the order of approximation we are using, by 



v 

T 



(aA) 2 
8tt 2 



1 



1 + 



A, 



16tt 2 



Ar 



(5.9) 



where Z^i = 1 — ^Ar/2 + O(Ap) is the multiplica- 
tive renormalization factor and c(t) = c — hi(a(t)m) . 
The renormalized coupling constant Ar is defined by 
A = Z x Xr, with Z x = 1 + §Ar/ 2 + 0(A|). With these 
definitions, we indeed obtain a renormalized equation of 
motion 



(d 2 + M 2 



vr. 



Ar/[ 

2 Vr 



o(x R ) 



3A R 7J 



0(A 



(5.10) 



and a renormalized gap equation 



1 



1 



M 2 {t) = -vi--rh 2 a\t) 1 + 



Ar 
16tt 2 



0(A r ) , 



(5.11) 
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where 

jf (t) = j£(t) - [(a(t)A) 2 - M 2 (t) In (o(t)A)l /(8ir 2 ) . 

(5.12) 

If our assumptions about the large-fc behaviour of the 
integrals are correct, then these are finite relations be- 
tween renormalized quantities. (To simplify matters, we 
have neglected a term in ( |5.10| ) involving the time deriva- 
tives of I2 which is finite, but probably an artifact of our 
regularization procedure.) 

For the symmetric state which exists prior to to, the 
renormalized mass fj,(t) is defined in a similar manner to 
M(t) and satifies the gap equation 



Ar f, 



^(t) = '-^lf{t)-~a\t)rh* 1 



Ar 
16tt 2 



- o{\l ) , 

(5.13) 

where lf(t) is defined in the same way as j£(t), but using 
the 4> propagator and mass n(t). 

C. Initial conditions at t = to 



The equation of motion ( |5.10| ) needs, of course, the 
initial values ur(£o) and Or (to) and these are obtained 
from the continuity of (<p 2 (t)) and d(4> 2 (t)) / dt at t = to. 
This expectation value is given by (4> 2 (t)) — if + 0(A) 
just before to and by definition it is equal to X~ 1 v 2 (t) 
just after to- In terms o|wr, the continuity conditions 
are (in the notation of fl4.8|)) 



(5.14) 



v 2 R (t) = X R lf(t) + 0(X 2 R ) 



As a consequence, we find that M(to) = /i(to) a t this 
order, suggesting the satisfactory possibility that the ex- 
citations of the state at to described by (f> and by £ are 
essentially the same. 

Our final task is to investigate how the initial con- 
ditions (HH) and (p~25l ) for the generalized occupation 
numbers Nu(t) and iVf(t) might be renormalized. It 
turns out that a far-reaching modification of the strat- 
egy u sed to obtain these conditions is needed. Indeed, 



( 4.23 ) has a potentially disastrous conseqence, since it 
implies that 



9 C n(to,to;k) = -f C k (t )ff(to) [N C k (t Q ) + JVf (t„) 



Ro(k) 
~ 2U 2 ■ 

On integrating this relation we discover that 



d 3 k 



9n (to, to; k) 



1 

2Ul 



d 3 k 
(27p 



2^T(*o). 



(5.15) 



R Q (k) 
(5.16) 



The renormalization programme we have outlined so far 
clearly requires that the ultraviolet divergences of i£ and 



If should be the same, which is incompatible with (5.16). 
The problem arises from the fact that we need to match 
two different perturbation series at to, whereas a con- 
sistent renormalization scheme can be implemented only 
within a single systematic expansion. In view of this, 



we replace <?£g(i, t') in (|4.19| ) and ( |4.2l| ) with a function 



g^a(t, t 1 ) obtained by summing those contributions to the 
2-point function for C, which are of order A , when v is 
regarded as of order A 1 / 2 , as it is in the </>-theory. This 
function is 



v(t) 



2v(t) 2V ' 'v{t) 



v(t') 



where 



5{t)= v 2 {t)-\ll(t) 



1/2 



(5.17) 



(5.18) 



and v{t) is understood to be evaluated at lo west order. 
On evaluating the modified version of (4.19) at t' = to, 
we find 



cr2 (*o)gii(to,to;fc) = ^ 



it (to, t ; k)- 4(to,t ;fc) 



(5.19) 



and we want to find initial conditions for for which 
this condition holds. It clearly will hold if a(to) — and 
if (to, to; k) = I2 (to, to; k), and this is in fact the only 
possibility. On integrating ( 5.1 9| ) and using in ( 5.18| ) the 
fact that v 2 (to) = XI f (to), we obtain 



7f(t )-7f(t ) 7f (to) = ^ ^(to) 2 -^^) 2 , (5.20) 



which implies that if (to) = Ii(to) and hence a(to) = 0. 
Now, we found above that the <p and £ masses /j,(t) and 
M(t) are equal (at lowest order) at t = to, so we can 
choose the same mode functi ons f k(t) for both propaga- 
tors near t = to- Evidently, ( 5.19| ) is satisfied if we also 
choose Ni,(tp ) — Nf (tp) and, in fact, the modified ver- 
sions of (f4.18| ) - ( 4^21 ) all hold at t = t' = t if we choose 



^fc(^o) = A^(to). We thus reach the apparently sat- 
isfactory conclusion that the propagators g^a(t, t') and 

5^(t,t') and their first derivatives are the same at to, 
so that the cf> and C excitations are essentially the same. 
The requirement ( [2.1£ ) is of course satisfied for N^,(t), 



since it is already satisfied for N^(t). The new bound- 
ary conditions do, of course, spoil the exact cancellation 
between the composite propagator of figu re 2(h) and the 
unwanted K.5 discussed in section [V B 1 . This propaga- 
t or ca n be calculated explicitly using the representation 
( |2.13 ) for the internal lines and setting the anchored time 
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arguments equal to to after the required derivatives have 
been performed. We find that derivatives of 9(t — t') pro- 
duce terms proportional to 5(0) which are again exactly 
cancelled by IC5, but there are no w resid ual terms. These 
contribute to the function B a p in (4.31) and do not enter 
into our explicit calculations. 

At this point, then, we have a fully renormalized equa- 
tion of motion and gap equation and finite initial condi- 
tions for both the equation of motion and the propagator 
g a g- To the first non-trivial order of our approximation 
scheme, all the ingredients are to hand for the calculation 
of the expectation values we want, which must necessar- 
ily be completed by numerical methods. This apparently 
satisfactory state of affairs conceals, however, an uncom- 
fortable feature of perturbation theory which must now 
be exposed. If we define 



ol(t)=vl(t)-\ K ll(t) + 0(\i) 



(5.21) 



then ( |5.1CD and ( PH) are consistent at this order with 
the pair of equations 



■& K + M 2 (t)a K = \al + 0(\l), 



(5.22) 



qualitatively in the expected manner. This ambiguity is 
inherent in any perturbative approach to our problem, 
which necessarily requires us to match the two different 
expansions about the states indicated schematically in 
figures 1(c) and 1(d). Intu itively, we would like to in- 
terpr et th e definition ( 5.21 ) as yielding or = (<fi) = 
and (5.24) as approximating the locations of the most 
probable values of (f> in figure 1(d). In practical terms, 
the most appropriate strategy seems to be that set out 
earlier. Namely, we make no reference to <tr, but deal 
exclusively with vr, which can be calculated both in the 
(^-theory and in the (-theory. Within the (-theory, we 
adhere to a strict series expansion of ( 5.1 0| ) and (5.11) 
in powers of Ar, treating the initial values obtained from 
the </>-theory as purely numerical input to these equa- 
tions, even though the numbers are obtained within the 
0-theory from an expansion in powers of Ar. In principle, 
at least, this yields a sequence of approximations which 
can be pursued systematically to arbitrarily high orders. 



VI. DISCUSSION 



M 2 (t) = \a 2 K (t) + iA R /f (t) 

-\m 2 a 2 (t) [1 + ±A R c(t)] + 0(\l) . (5.23) 

Suppose, for the sake of argument, that a(t) approaches 
the constant value a — 1 a t lat e tim es. There are 
two steady-state solutions to ( 5.22 ) and ( 5.23| ) to which 
the final state might correspond. The first is or = 
V3M + 0(A|) and M 2 = m 2 - \\ K l{ 



0(X\). 



This 



corresponds as in section [I B to a superposition of states 
in which ((f)) = ±or, and is the state that we might hope 
to see emerging. The second solution is or = and 



M 2 = -\m 2 + |A R /f + 0(A R ). In terms of the con- 



ventional description of spontaneous symmetry breaking 
in Minkowski spacetime, this of course corresponds to 
the ill-defined perturbation theory with ((f)) — 0. Which 
of these situations will emerge from our dynamical de- 
scription depends crucially on h ow we treat our pertur- 
bative expansions. Using ( [5.21 ) as it stands, the initial 
conditions we have derived would imply that or = 
for all times (or being essenti ally a renormalized ver- 
sion of the a introduced in fl5.1§| )), and hence that 
AR-ff at all times. Although we have 



Ar/ : ' 



only first-order results in hand, this strongly suggests 
that that the (-theory would then simply reproduce the 
</>-theory, but by a more complicated route. This pertur- 
bation theory is unsatisfactory when wr becomes of order 
1, as it eventually must do. On t he other hand, if we use 
the actual O(Ar) result in ( 5.1C ) to define 



cr = v R 



(5.24) 



and truncate the expansion at this order, then or and 
or receive nonzero initial values at to and should evolve 



In this paper, we have proposed a perturbative approx- 
imation scheme through which quantum-field-theoretic 
expectation values might be estimated during the course 
of a symmetry-breaking phase transition. The scheme 
has two distinctive features, both of which seem to be 
essential to this problem. First, in order that low-order 
calculations adequately reflect the evolving nonequilib- 
rium state, this state is described in terms of its own 
quasiparticle excitations. Lowest-order propagators de- 
scribe the propagation of these excitations and incorpo- 
rate approximations to their mass (or, more generally, 
their dispersion relation) and damping rate, which can be 
obtained as the solutions of appropriate self-consistency 
conditions. Second, because we deal with an exact sym- 
metry 4> <-> —(f), which is a symmetry both of the Hamilto- 
nian that governs time evolution and of the initial state, 
the scalar field <f> can never acquire a nonzero expectation 
value. In fact, the symmetry can never, properly speak- 
ing, be broken at all. In Minkowski spacetime, neverthe- 
less, the effects conventionally associated with a nonzero 
value of (<f)) can be recovered in perturbation theory by 
dealing instead with ((f) 2 ) and we have here generalised 
this description to accommodate the state emerging dy- 
namically from the phase transition. (As noted in section 
U the value of ((f)) does not in itself provide a full char- 
acterisation of the noncquilibrium state. Nor, of course, 
does the value of ((f) 2 ), but this quantity plays an in- 
dispensable role in the perturbative description of this 
state.) 

Inherent in any attempt to apply perturbation theory 
to this problem is the difficulty that ((f) 2 ) is formally of 
order A in the early, symmetrical state, but of order A -1 
in the state that emerges from the transition. This nec- 
essarily means that two inequivalent perturbation series 
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must be used, before and after a time that we have called 
to- Intuitively, we suppose that before to, the most prob- 
able value of <j> is zero, while after to the most probable 
values are, say, ±cr/vA. However, perturbation theory 
gives no clear indication of the instant at which this bi- 
furcation occurs, so in practice an ad hoc choice of to, 
perhaps optimising in some way the matching of the two 
approximations, would be required. More importantly, 
perhaps, the above difficulty leads to two specific short- 
comings of our analysis which deserve some emphasis. 
One is that we had to perform in section III B a somewhat 



singular transformation of the path integral from which 
expectation values are computed, which led to ambigu- 
ities concerning the treatment of total time derivatives 
and the meaning of derivatives of propagators when their 
time arguments coincide at to ■ We succeeded in resolving 
these ambiguities sufficiently for the immediate purposes 
of the explicit calculations described in sections IV and 
0, but have not obtained any well-founded general reso- 
lution. It is possible that a sufficiently careful derivation 
might avoid these ambiguities, but we have not yet been 
able to achieve this. The second shortcoming is that 



discussed in section V C , namely that our perturbative 



estimate of Vn(t) fro m the equation of motion (5.1C) and 



the initial contitions (5.14) depend crucially on how the 
perturbation series is organised. We gave a prescription 
which should apparently produce the expected type of 
behaviour, and can be pursued systematically, but it is 
not unique. We do not think that any significantly better 
strategy is available within perturbation theory. Possibly, 
some non-perturbative means of improving the matching 
between the two expansions might be found, but we have 
not yet been able to do this. 

Quite extensive numerical calculations are needed to 
extract concrete results from the perturbation theory we 
have constructed, and we hope to report such results in 
future publications. In the absence of further approxi- 
mations beyond the perturbation expansion, the numer- 
ical problem involves the simultaneous solution of sev- 
eral coupl ed nonlinear equati ons: the equation o f mo - 
tion ( 5.10 ); the gap equation ( 5.11 ); the equations ( 2.15| ) 
and (2.17) for mode functions and generalized occupation 
numbers, and two equations we have not given explicitly 
for the dissipative coefficients ak(t) and 7fe(t), which in- 
volve at least 2- loop integrals, as detailed in |2l]]. If the 
scale factor a(t) is also to be determined self-consistently, 
then there are also two independent field equations con- 
tained in (1.1). The nature of the solution is rather dif- 
ficult to forecast, but we wish to speculate briefly on 
some particular features, partly in the light of investi- 
gations reported by Boyanovsky et. al. for the some- 
what simpler case of the large-iV model. These authors 
studied a phase transition roughly of the kind envisaged 
here in |31[ | using a fixed de Sitter background and in 
JlT| by solving the Friedmann equation simultaneously 
with the field theory problem. The field-theory equa- 
tions th ey so lve ha ve, to some extent, a similar struc- 
ture to ( 5.22| ) and ( |5.23| ). There are trivial differences 



arisi ng fr om the fact that, in order to obtain the ac- 
tion ( |l.2|) (which is equivalent to a Minkowski-space the- 
ory with a time-dependent mass simply proportional to 
a(t)) we have used <fr(x) — a(t)Q(x), where is a 

genuine scalar field under general coordinate transfor- 
mations, taken this field to be conformally (rather than 
minimally) coupled to gravity and used a conformal time 
coordinate t. There are, however, more fundamental dif- 
ferences. Most obviously, the one-loop equations are ex- 
act for the large- N model, whereas those we have given 
explicitly are just the lowest non-trivial order of our per- 
turbative approximation scheme. In the large-iV model, 
the quantity denoted in [ pdpl| ] by rj(t), which is analo- 
gous to our aji(t) is the expectation value of one of the 
N fields. The evolution of this field is entirely classical, 
while the quantum modes contributing to the self-energy 
integral I\ arise purely from the N — 1 transverse fields. 
As a consequence, the mode equation corresponding to 
( |5.l|) is identical for k = to the equa tion of motion for 
i](t), while here it coincides with (5.22) only when aji(t) 
is small. In the absence of explicit symmetry breaking, 
one has rj(t) =0 and the calculation is analogous to using 
our 0-theory for all times. The large- N limit produces, of 
course, a constrained Gaussian model, in which there is 
no s cattering, so the dissipative effects treated in section 
II A are entirely absent. 

Suppose, then, that our 0-theory were u sed f or all 
times. The relevant gap equation is then (|5.13|) and 
the mode functions contained in lf(t) obey ( p!l| ) with 
M(t) replaced by /i(t). Initially, is positive and if (t) 
should vary only slowly with time. Since a(t) increases 
with time, fj, 2 becomes negative, at which point the mode 
functions fk(t) for which k 2 + /i 2 < begin to grow, and 
so therefore does lf(t). If this growth is fast enough, 
then /i 2 might again become positive and, perhaps, a fi- 
nal state might emerge in which p 2 remains positive. In 
the large- TV model, this is more or less what happens, ex- 
cept that the quantity analogous to /i 2 approaches a van- 
ishingly small value, corresponding to a universe popu- 
lated entirely by Goldstone bosons. In our N — 1 model, 
a late-time state with fi 2 > would make little sense 
(though this does not rule it out as a solution to our 1- 
loop equations!). In such a state, the thermal contribu- 
tion to the effective mass is still large enough for the sym- 
metry to remain unbroken (in the conventional sense): 
in the (not strictly appropriate) language of equilibrium 
field theory, the system would have reheated to above its 
critical temperature. In order to obtain a sensible final 
state, it is essential to move to the C description at some 
point durin g the growth of the unstable modes. From 
( |5.22 ) and ( 5.23| ), it is clear that when or is positive, it 
grows more rapidly than any of the unstable modes (and 
the same should be true of i>r in our preferred approx- 
imation scheme), so that M 2 is more likely than /i 2 to 
become positive at late times, as it should. 

Finally, we comment briefly on two significant issues 
raised by the growth of unstable modes. First, these 
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modes give rise to large propagators which threaten to 
make perturbation theory useless. However, as indicated 
above, the growth does not continue indefinitely, and the 
period of growth will become shorter if the coupling is 
made stronger. In fact, the results presented in |3l]] show 
a quantity analogous to XIi(t) growing to a value of 1 
in the large- TV model and to about 1/3 in a Hartree 
approximation which is more similar to our 1-loop ap- 
proximation. While these numbers are not small, they 
seem to offer hope that perturbation theory may not be- 
come grossly unreliable. In our ^-theory, moreover, ver- 
tices involving £ n carry factors proportional to v 2 ~ n (see 



( |3.50[) ). Since each factor of ( in £? nt corresponds roughly 



to one mode function fk(t) inside a Feynman diagram, 
the growth of these mode functions will be largely off- 
set by the growth of i>r. The second issue concerns the 
nature of the state at late times. In jy], the small- fc 
modes of the large- N model are interpreted as constitut- 
ing an effective classical field <j) c g ~ (0 2 )r. Here, this 
role is explicitly assumed by vr, but the remaining quan- 
tum field £ still contains the growing modes. Once these 
mode s h ave acquired large amplitudes, the mode equa- 
tion (5.1) provides no mechanism for these amplitudes to 
decay, and at first sight, this seems at odds with the per- 
turbative description of the thermalized state one might 
expect finally to emerge. One would expect to describe 
this state in terms of mode functions with amplitudes 
of order (fc 2 + M 2 ) -1 / 4 , and the presence of large am- 
plitude modes seems to suggest coherent behaviour that 
ought to have been absorbed in the expectation value 
ur. However, the decomposition of our propagators g a p 
into mode functions and generalized occupation numbers 
iVfe(t) is not unique. At late times, it will be possible to 
reexpress these propagators in terms of small-amplitude 
mode functions and large occupation numbers for the 
previously unstable modes. We then see that the dissi- 



pative mechanisms of section II A should cause the ex- 
pected thermalization (though this does not depend on 
how we choose to write the propagators down!). Indeed, 
dissipation will be present during the period of instabil- 
ity as well. Its effects are difficult to forecast, but some 
inhibition of the growth of unstable fluctuations seems 
likely. 
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FIG. 1. Schematic representation of a time dependent ef- 
fective potential and of the probability density for (f>. 

FIG. 2. Diagrammatic representations of the propagators. 
Solid lines (a) represent g^, dashed lines (b) represent g* b 
and dotted lines are cj> propagators whose time arguments are 
both equal to to- 

FIG. 3. Vertices representing interactions in the <f> and ( 
theories. 

FIG. 4. Low-order diagrams contributing to (C(t)). An- 
chored contributions (a) contain composite propagators with 
vertices fixed at to while the free contributions (b) contain 
none of these vertices. 

FIG. 5. Leading-order contributions to {cj> 2 (t)4> 2 (t')) when 
(a) both times are earlier than to, (b) one time is later than 
to and (c) both times are later than to- 

FIG. 6. Low-order contributions to the connected 
two-point function {Ga(t)£p(t')}. 
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